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Abstract. For more than a century Liesegang patterns - self-organized, quasi- 
periodic structures occurring in diffusion-limited chemical reactions with two com- 
ponents - have been attracting scientists. The pattern formation can be described 
by four basic empirical laws. In addition to many experiments, several models 
have been devised to understand the formation of the bands and rings. Here we 
review the most important models and complement them with detailed three- 
dimensional lattice-gas simulations. We show how the mean-field predictions can 
be reconciled with experimental data by a redefinition of the distances suggested 
by our lattice-gas simulations. 

1 Introduction 

In recent years the general interest in self-organized structures is growing, triggered by the 
idea of cheap and fast production of nano-scaled devices. One of the promising effects for 
obtaining such devices is Liesegang pattern formation, based on a reaction-diffusion process. 
Experimental evidences for Liesegang patterns in solid materials on the nano-scale have already 
been obtained, see, e. g., |l|2|3j . For example, periodic patterns of silver nano particles in glass 
were observed by electron microscopy pQ. If it became possible to control the growth of such 
patterns with experimentally tunable parameters, self-organized optical devices could be made. 
In addition, the Liesegang phenomenon is an interesting research topic on its own due to the 
simple patterns arising out of complicated reaction and diffusion processes. 

Since the first description of Liesegang rings in gels by the German chemist Raphael Eduard 
Liesegang in 1896 [4 , many experiments and several models have been devised to understand 
the formation of the bands and rings. Although the basic problems are solved and four universal 
empirical laws have been found common to all experimentally observed Liesegang phenomena, 
there are still open questions. In literature misunderstandings regarding the exact definition of 
the measured quantities cause some discrepancies between the experimental results and the pre- 
dictions of theoretical (mean-field) models. In this paper, we review the most important models 
for Liesegang pattern formation and complement them with extensive lattice-gas simulations 
of the reaction-diffusion processes. Based on our three-dimensional (3d) simulation results, we 
find a way to reconcile most experimental observations with the results obtained in mean-field 
models. In addition, we obtain evidence that fluctuations in 3d reaction-diffusion processes may 
play an important role in Liesegang pattern formation on the nano-scale. 

The paper is structured as follows. Section [2] describes the main experimental findings, 
which can be summed up in four empirical universal laws describing the Liesegang patterns. 
Section [3] reviews previous work on models reproducing these laws. After an overview including 
a discussion of recent trends, we focus first on reaction-diffusion models with thresholds, and 
describe afterwards the spinodal decomposition model and a kinetic Ising model. Section [4] is 
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Fig. 1. (colour online) Experimental examples of Liesegang patterns in test tubes. Agent A is injected 
from the open end (top) of the tube, which contains agent B dissolved in a gel, yielding a linear geometry 
(Figure by D. B. Siano, http://commons.wikimedia.0rg/wiki/Image:Liesegangrings.jpg). 

devoted to our studies based on lattice-gas simulations, presenting both the numerical method, 
the main findings, and a suggestion for reconciling the mean-field predictions with both, the 
experimental data and our 3d simulation results by a redefinition of the distances. Section [5] 
summarizes our findings and gives an outlook on further work on Liesegang pattern formation. 



2 Experimental findings and empirical universal laws 

In his original experiment, Liesegang covered a glass plate with a layer of gelatin impregnated 
with potassium chromate [4]. Then he added a small drop of silver nitrate in the centre. As a 
result, silver chromate was precipitated in the form of a series of concentric rings with regularly 
varying spacings. These rings became famous as Liesegang rings or more generally Liesegang 
patterns. Shortly after the first experiments, Wilhelm Ostwald presented an explanation for 
the occurrence of the rings [5], which is still the basis for most of the models today. The 
next important experimental findings followed several years later. In 1903, Morse and Pierce [6] 
investigated the formation time of the bands, observing diffusional dynamics, i. e., the time law. 
Jablczynski [7] showed twenty years later that Liesegang patterns follow a geometric series, i. e., 
the spacing law. Based on this observation Matalon and Packter |8|9j investigated the functional 
dependence of the positions of the bands on the concentrations of the reacting agents in 1955. 

The geometry of the pattern depends on the initial conditions for the reacting agents. Usually 
one agent (the inner electrolyte), represented by the B particles in the models, is initially 
homogeneously distributed in the sample or gel. The second agent (the outer electrolyte), 
represented by the A particles, is injected. If A is injected in the centre, precipitation rings are 
formed. If B is homogeneously distributed in a cylindrical tube and A is injected from one end 
of the tube, bands form perpendicular to the motion of the reaction front, see Fig. |TJ The second 
experimental setup is more appropriate for a theoretical description because it is effectively one 
dimensional (Id) [10]. Most of the theoretical and experimental work was done using this linear 
configuration; here we also focus on the 'band' setup. Although Liesegang rings are basically a 
projection of the bands onto polar coordinates, spiral patterns have been observed in the ring 
configuration exclusively. They have no equivalent in the linear configuration. 

To define the basic observables, we assume that the nth band forms at time t n at distance 
x n from the side where A is injected. The width of the nth band is denoted by w n . Firstly, the 
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Fig. 2. Results of a 3d lattice-gas simulation of Liesegang band formation, (a) Injecting A from the 
left, Liesegang bands (black) with increasing distance form in the simulated lattice of 32 x 32 x 2048 
sites. To obtain the projection, all 32 two-dimensional slices were placed next to each other vertically, 
(b) The ratio x n+ i/x n of the positions of the (n + l)st and nth band is plotted versus n. See Fig. 6 in 
Sect. 4.3 for the parameters of the simulation. 



position x n is empirically found to be proportional to the square-root of the time t n 

t >. 



(1) 



This rule called time law in literature has been confirmed experimentally many times |6|11|12|13|14I15|16I17|18| . 

Figure [5] shows that the position x n of the nth band follows approximately a geometric 
series converging to the spacing law [7] 



x n oc (1 +p) n <^ x n+ i/x n —>■ 1 + p for large n 



(2) 



with p > the spacing factor. Typical empirical values for p reported in literature range from 
0.05 to 0.4 [112113116118119] . 

A first systematic experimental analysis of the functional dependence of p on the concentra- 
tions of agents A and B was done by Matalon and Packter [819] . They gathered experimental 
results on p and found a functional dependence on the concentrations of both agents, do and 
&o, respectively. This dependency is known as the Matalon- Packter law and takes the form 



p = F(b ) + G(b )-^ 
a Q 



(3) 



with F and G as dimcnsionless monotonously decreasing functions of bo [20) ■ F and G are not 
dependent on ao making p linearly dependent on 1/ao for constant bo. This can serve as a test 
for the validity of the Matalon-Packter law. In Sect. 13.3.31 we will show that Eq. |3]) is just 
an approximation of a more general law in the limit of reaction fronts much faster than the 
diffusion of the B particles. This limit holds for most of the experiments, because ao is usually 
much larger then bo while both electrolytes diffuse equally fast. An exception seems to be the 
experiment with nanoscale silver particles in glass [1] where the diffusion of silver ions, the inner 
electrolyte, might be one magnitude faster than the diffusion of the outer electrolyte, hydrogen. 
This leads to a slower motion of the reaction front compared to the diffusion of the silver ions. 

As can be seen in Fig.[2]Ja) the Liesegang bands are getting broader for larger n. It is possible 
to set up an empirical law describing the width w n of the nth band as function of the position 
x n . Experimentalists use two competing versions of the width law, 



w n = p,ix n + u 2 



and 



w n cx x. 



(4) 
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see, e. g., |12|19|21] and [IB 22 23 24J, respectively. We are not aware of any papers comparing 
the two versions. In addition, different values of the exponent a in the second version have 
been published. The experiment with nanoscale silver particles in glass can be fitted by a ~ 0.7 
P] (using the second version of Eq. J4}), but the results can be equally well fitted by the first 
version. Most other papers report larger values of a, e. g., 0.9 < a < 1 [Mj. Publications of 
early reactive lattice-gas simulations, where the second version was first introduced, fit best with 
values of a ~ 0.5 — 0.6 22.25], but no comparison with the first version of Eq. (j4j was done. It 
is difficult to distinguish between both versions of the width law because the number of bands 
is limited and the widths w n of the bands are small compared to their positions x n . Theoretical 
works prefer the second version with a = 1 [18 23 24]. In Sect. 14.51 the different versions will be 
tested on our results of lattice-gas simulations. We will also propose an alternative, theoretically 
well grounded approach unifying both versions there. 

The time law (Jl]) is a simple consequence of a diffusion process. The position Xf(t) of the 
reaction front between A and B moves proportional to \fi with a prefactor that depends on ao 
and bo [IS]. The other three laws cannot be explained so easily; one needs models describing 
nucleation and growth of the bands. These models can be categorized into two types. The 
first type is based on diffusion and reaction dynamics plus some thresholds which account for 
nucleation and growth. There exists a wide range of modifications. The second type of models 
uses well established phase separation techniques to explain the pattern formation based on 
spinodal decompositions [57] or a kinetic Ising model [38] • All of these models will be reviewed 
in the next section. 



3 Models 

3.1 Overview and recent trends 

Reaction-diffusion models with thresholds are the oldest models describing Liesegang pattern 
formation. Only a few month after the first experiments by Liesegang, Ostwald gave a pos- 
sible explanation of the pattern formation process on the basis of supersaturating liquids 5 . 
He suggested that the precipitation is not a result of a balanced reaction but must happen 
spontaneously when the concentration product K of the reactive partners reaches a critical 
concentration K sp . The precipitation grows until K falls under a second concentration K p . 
Both thresholds set the stage for several models differing in the detailed description of the 
reaction process. In principle they all follow a reaction scenario 

A + B^---C >D, (5) 

and differ in the way the intermediate stage ■ • ■ C ■ ■ ■ and the precipitation — > D are described. 
For the concentrations of the different chemical agents A, B, . . . the symbols a,b, . . . will be 
used in the following. 

The empirical laws described in the previous section became the basis for the theoretical 
work which started in the fifties. First analytical descriptions |29|30|31j showed that the pat- 
terns can be explained by diffusion and reaction processes with a moving reaction front. When 
numerical simulation became feasible about 90 years after the first experiments it was possible 
to reproduce patterns dynamically by studying the reaction-diffusion equations [14121] . The 
fact that a novel approach was introduced in 1999 [23127] shows that the Liesegang pattern 
phenomena are still an active research field in both, experiment and theory. 

Apart from the novel simulation approach the trends in the literature follow two general 
lines. The first line works on an analytical description of the basic laws [24 32 33 or even looks 
for more general laws [34] . The second line varies the initial conditions and the geometrical con- 
figurations of the setup to understand the phenomena better and to find ways for applying the 
pattern formation in interesting engineering problems. Alternative geometries [35] and complex 
3d boundary conditions [2 3 36J have been studied experimentally. Variation of the reaction 
terms [37] , additional terms for a dissolution of the bands [15138] , as well as open systems [39] 
and, last but not least, systems with an additional electric field [15|40|41(42] have been studied 
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to vary the patterns. In very recent work the motion of the reaction front is even detached 
from the phase separation, which might be initiated, e. g., by a temperature gradient |43j . One 
hopes to understand and control the pattern formation process such that bands with constant 
distance can be generated, see |2|3j for experimental work in this direction. 

Most of the theoretical work is based on an analytical study or numerical solution of dif- 
ferential equations. Since such equations are always based on the concentrations a, b, . . . of the 
reacting agents, they generally yield mean-field solutions. Although such mean-field solutions 
can reproduce the basic laws listed in the previous section they have two disadvantages. Firstly, 
they cannot account for the statistical character of the reactions. Hence, the influence of particle 
number fluctuations and thus the stability of the patterns cannot be investigated adequately. 
Although there are some ideas to include fluctuation by an additional noise term in mean-field 
models |2|40j . we think that models with intrinsically statistical character are more adequate. 
This is particularly true for mesoscopic and nanoscopic systems, where fluctuations become 
more important. First experiments on microscopic scales indicate that fluctuations might play 
an important role [l|2|3j . Secondly, the differential equations may not describe the microscopi- 
cal structure of the bands. Chopard et al. 22 25J proposed an alternative approach by taking 
the basic principles of the mean-field description and implementing them in a reactive lattice- 
gas simulation. Such simulations can be very helpful in reconciling mean-field predictions with 
experimental data, as we will see in Sect.0] A similar approach seems also possible using Ising 
models [28] , 

In the following sections, the mentioned quantitative models and simulations will be re- 
viewed, except for the lattice-gas simulations to be presented in Sect. HJ 

3.2 Ion product saturation models 

In the first and easiest quantitative models for Liesegang pattern formation, the ion product 
models, the precipitation takes place without an intermediate stage; i. e., there is no C stage in 
© |29I30I31| . Like in the Ostwald model, nucleation A + B — > D occurs when the local product 
concentration K — ab reaches the threshold K sp . The precipitates of D will grow and deplete 
their surrounding of A and B. As the reaction front proceeds, the product concentration around 
the immobile precipitate decays until growth becomes impossible. Wagner [29] , Prager [30] and 
Zeldovich et al. [31] could show that these ingredients yield patterns which obey the time law 
(P) and spacing law The diffusion profiles of a and b are described by a system of coupled 
integro-differential equations with given boundary conditions. These equations can be rewritten 
in a more convenient way, if local coordinates are introduced in the form A„ = x n /(2y/Dt n ) 
and 7„ = x n /x n -\. Taking n — > oo and xq — > renders the equations mathematically solvable. 
The local coordinates A„ — > A and j n — ► 7 yield the time law and the spacing law, respectively. 
However, since these laws were the result of a continuous limit, the dynamics of the process was 
lost. Furthermore, the formation of the Liesegang bands was taken for granted and not proved 
by solving reaction-diffusion equations. 

The first dynamical version of the ion product model proposed by Ross et al. [14] was 
deduced from the elementary chemical reaction A + B — > D and solved numerically. Besides 
the typical diffusion terms for A and B, 



with diffusivities Da and Dg, respectively, the equations include a reaction rate R constructed 
for a specific experimental configuration. In principle, R mimics an auto-catalytic growth pro- 
cess with a growth rate proportional to an increasing function of the supersaturation S = 
(ab/ K sp ) — 1. The authors tested different power laws and exponential functions and observed 
pattern formation for a very quickly growing function R(S) only. Beside the confirmation of 
the time law no quantitative observations of the other laws were reported. 

Later Zrinyi et al. [44] pointed out that no detailed description of the agent transport is 
needed, since precipitation takes place on a faster time scale. It is thus possible to model the 



dta = DAd^.a — R and 
d t b = D B d 2 x b-R 



(6) 
(7) 
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growth process by a given critical threshold K p . The corresponding source term for additional 
D can be written as a Heaviside step function, c\(a, b, d) 0(ab — K p ). The growth starts when 
the product concentration exceeds a second threshold K sp > K p , where the rate is independent 
of d. The partial differential equation for the precipitate concentration d can thus be written 

as pres] 

R = d t d = a(a, b, d) 0(ab - K p ) + c 2 0(ab - K sp ), (8) 

where c\ and c 2 are model-dependent functions of the parameters (initial concentrations and 
diffusion coefficients), and c\ also depends on a, b and d. Liesegang patterns observed in a 
wide range of systems can be modelled this way. In particular, it was possible to show that the 
time law and the spacing law hold for a large set of parameters [44] . The Matalon-Packter law 
was confirmed only if the initial concentration of the outer electrolyte clq is sufficiently large. 
However, p was shown to stay a monotonously decreasing function of ao even if ao is small. 
Furthermore, p increases monotonously with K sp in a non-linear way, and it is anti-proportional 
to the diffusivity Da of the outer electrolyte. The results for a specific p do not depend on the 
individual parameters but on ratios of parameters with equal dimension - a simple consequence 
of the mean-field character of Eqs. © to ([5]). The findings are consistent with an experimental 
categorisation of the patterns by the product and the difference of the initial concentrations 
[13]. The width law, however, was not investigated. 

In closely related approaches, growth was modelled proportional to the concentration of 
the precipitate, corresponding to K p — and c\ oc d in Eq. ([8]) 32 33J . Focusing on the 
Matalon-Packter law ([3]), Antal et al. found that F is constant and G oc y/ Db j Da K sp /b^ 
|32j . This result was achieved analytically based on the assumption bo <C ao; hence the failure 
of the Matalon-Packter law for low ao could not be observed. This result is inconsistent with 
p ex I/Da in [44], due to different assumptions for c\ and c 2 . Lebedeva et al. searched for a 
criterion for the pattern formation [33], finding <P = K sp k 2 D b / 'u 2 (t) > 2 + with k^ the 
growth rate and u(t) the velocity of the front (decreasing with time). Patterns are thus formed 
when the nucleation threshold K sp is high and the growth rate k 2 is large. The time dependence 
of u{t) explains why patterns do not start at xo — but an empty ("plug") zone exists at the 
edge. The first two results match with those of Ross et al. [13] and Zrinyi et al. [44] . who did 
not study the velocity of the reaction front. 

3.3 Nucleation and growth models 

The next models with slightly increased complexity separate the reaction-diffusion process 
A + B — > C from the nucleation and growth process C — > D, introducing an intermediate 
state C, see Eq. ([5]). This, however, considerably simplifies the analysis. In the literature such 
models are often called nucleation and growth models because the formation of D depends on 
the nucleation and accumulation of C rather than on the product concentration of A and B. 

This model was first introduced by Keller et al. |46] who analysed reaction-diffusion equa- 
tions and confirmed the time law ([1]) and the spacing law ([2]). In contrast to Wagner and Prager 
[29|30j they could calculate the positions x n of the bands without a priori assuming band for- 
mation and stopping of the band growth. The first numerical solutions for a nucleation and 
growth model were presented by Dee [21]. He used a similar technique as Ross et al. [14] (see 
previous section) and determined the nucleation and growth criterion by classical nucleation 
theory A detailed work concerning the nucleation and growth of silver particles similar to 
Dee was presented recently [47] . Since the reaction-diffusion process can be separated from the 
precipitation process, we will discuss the two stages separately. 

3.3.1 First process: reaction-diffusion A + B — > C 

The simple reaction-diffusion process A + B — > C is important independently of Liesegang 
patterns since it is a basic process in many chemical reactions. To apply the results to the 
Liesegang pattern phenomena we will focus on a quasi- Id geometry where the reacting particles 
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Fig. 3. (colour online) Illustration of the concentration profiles obtained in the reaction-diffusion 
process A + B — » C without (a,b) and with (c) precipitation C — » D. The concentrations a{x,t) 
(red) and b(x,t) (green) are shown for fixed time t together with (a) the reaction rate R(x,i) oc ab 
(black) following Eq. Q, (b) the accumulated concentration c{x,t) (magenta), and (c) the rescaled 
concentration of the precipitate D (blue). Note the constant value of c(x,t) — Co behind the reaction 
front in (b) according to Eq. The dashed line in (c) indicates the threshold K sp ; the next band is 
started when c(x,t) reaches K sp . 



are separated at time to — with initial concentration ao and &o- The concentration profiles 
a(x,t) and b(x,t) for x > and t > are determined by the reaction-diffusion equations ([6]) 
and ([7]) with the initial conditions a{x, 0) = ao<9(— x), b(x,0) = boO(x) The reaction term 
R is again assumed to be proportional to the concentration product K = ab. Under these 
assumptions it is possible to calculate a(x, t) and b(x, t) as well as the reaction rate R(x, t) oc ab 
asymptotically for large times t [25], 



R(x,t) ~A f S R 



Wf(t) 



(9) 



where Sr is approximately a Gaussian [27150] . The result is illustrated in Fig. [3Ja). R(x,t) 
describes the reaction front and reaches its maximum value at Xf which scales as Xf(t) ~ \ft. 
The width of the front scales as Wf(t) ~ V with 7 = 1/6. The production rate Aj of C 
at x = Xf is proportional to t^^ with £ = 2/3. A scaling analysis regarding the impact of 
fluctuations [51] leads to the conclusion that the critical dimension, above which fluctuations 
become unimportant, is d c = 2. In 2d logarithmic corrections apply and in Id the whole 
dynamics change dramatically [51j . 

It is possible to approximate the effective diffusion coefficient D/ for the reaction front, i. e., 
the prefactor in Xf — \/2D ft as well as the concentration cq of C behind the reaction front. This 
requires a few assumptions [55] which will be usually fulfilled in typical experimental setups. 
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The reaction front R(x,t) has to be confined in a region (reaction zone) Xf(t) — Wf(t)/2 < 
x < Xf(t) + Wf(t)/2 for all times t, disregarding all reactions outside this region. Then the 
concentration profile of each agent can be written as a solution of a diffusion equation 26J. 
Furthermore, Eqs. ([6]) and (|7|) must be approximated by a quasi-stationary solution in the 
region — (D^) 1 / 2 <C x <C (-Dsi) 1 / 2 . Under these conditions the effective diffusion coefficient 
Df is given by an analytic implicit equation [52]. 

with H(x) = [1 — erf(a;)] exp(x 2 ). Here, erf(x) is the Gaussian error function defined as erf(x) = 
— 7^ Jq exp (— z 2 ) dz. Evidently the velocity of the reaction front does not depend on each of 

the initial concentrations, but on their ratio ao/bo, and it scales in the same way as Da and 
Db, since no characteristic time and length scales exist. The same approximation also yields 
that the concentration of C behind the reaction front is constant (see Fig.[3]Jb)), given by [32] 

The interesting point of these results, for the Liesegang pattern phenomena, is that only Co 
and Df are essential parameters. Since both depend only on the initial concentrations and the 
diffusion coefficients, the band formation is independent of the reaction rate which is hard to 
measure. This does not imply however, that the width of the bands or the specific concentration 
profiles are independent of the reaction rate. 



3.3.2 Second process: precipitation C — > D 

Describing the A + B — > C process exclusively by Cq, Df, and the Gaussian reaction front f9}, 
we now focus on the precipitation step. Like Ross et al. [14] (see Sect. I3.2[) . Dee [21] modelled 
the precipitation process using the classical droplet theory. Droz et al. [24 showed that Dee's 
approach can be simplified into 

d t d(x,t) =ciN[c(x,t),d(x,t)]+c 2 9[c(x,t)- K sp ] with N{c,d) oc d(x,t), (12) 

corresponding to Eq. with K p = and c\ and ci constants or proportional to c. A typical 
result of such a simulation is shwon in Fig.^c). Using this model, Dee confirmed the spacing 
law ([2]) and found the linear form of the width law (|4|) . Due to limited computational resources 
he was restricted to one example with only six bands and obtained no information on the 
Matalon-Packter law ([3]). 

In Eq. (fi"2|) the growth of the precipitates is restricted to the points of the nucleated particles, 
see also Fig. [3jc) . To obtain bands with a macroscopic width w n one would need nucleation at 
every point in w n since no non-local term exist. Although no macroscopic width is reproduced 
in this model one sees that the hight of the bands, which can be interpreted as their mass 
or integrated particle number, seems to grow linearly. This is in agreement with a theoretical 
prediction which will be introduced in Sect. 14.51 

Alternatively it would be possible to allow non-local growth such that C can become D in 
a broader surrounding of an initial nucleation centre. A popular version of the non-local term 
is [2153] 

nX-\-dx 

N{c,d) = I d(x',t)dx' 

x — dx 



0[c{x,t)-K p }. (13) 



In this case, a second threshold K p is needed to stop the growth. This threshold is not needed 
in Eq. (fT2|) because local growth stops if nucleation stops. 
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3.3.3 Quantitative mean-field predictions for p and the Matalon-Packter law 

It is possible to calculate the functional dependence of p in the nucleation and growth model 
described above [32] . 

Do r - 



co _ 1 _ D c 



K sp 2D f 



(14) 



This surprisingly simple analytical solution determines p from the basic parameters D f and cq 
(from A + B — > C), Dc (for the diffusion of C), and the nucleation threshold K sp . As expected 
for a mean-field solution, Eq. I|14p depends only on dimensionless ratios of the parameters. The 
first term characterizes the velocity of C versus the velocity of the reaction front, representing 
a measure of how fast cq could be reached. If Df <C Dc, G will quickly leave the reaction front 
and diffuse away. In this case it takes longer to reach the desired threshold K sp , and p becomes 
larger. The second term makes sure that no precipitation occurs if K sp > cq. If K sp is close to 
Co it takes longer to reach K sp and p is also large. 

The structure of Eq. (fT4"]) will be similar if nucleation and growth is modelled somewhat 
differently. All mean-field solutions must scale, and patterns emerge only if K sp < cq. For this 
reason Eq. (fT4|) can serve as reference to test further properties of Liesegang pattern formation. 
One of these properties is the Matalon-Packter law, which is clearly inconsistent with Eq. (fT4")l . 
since it does not scale linearly with 1/ao- However, Antal et al. [32] could show that the Matalon- 
Packter law is a special case of Eq. (fl4|) . if 6o <C do- Then Db <C Df and the Gaussian error 
function can be approximated by erf (a;) « 1 — exp(— x 2 ) (l — — . . .) j \pKX [54] . Applying 
this to Eq. (JTT]) we obtain c «6 (1+ § a ) = b (1 + 757757)- In the range 0.05 < b /a < 0.1 
interesting for experiments, Antal et al. showed numerically that Da/Dj is linear in ao/bo- 
This leads to Da/Dj w v\ + z^&o/ao, yielding an approximation for cq, 

cq = b Q 1 + i^i — + i/ 2 - , (15) 



Da D a a 

with 1 + v\jf- and ^2-§j numbers of order one. Since bo/a -C 1 it is possible to neglect the 
last term in Eq. fl"5"]) . It is not possible to argue this way already in the linear approximation 
of Da/Dj because the ^s do not have the same magnitude. Inserting the approximation back 
into Eq. (| 14[) we obtain 

D c K sp v x D c K sp v 2 bo b b 

= F(b ) + G(b )—, (16) 



DA((Tibo — K sp ) D A {(Ti - K sp ) a a 

which is in the form of Eq. ([3]). The assumptions needed to derive Eq. (fTB)) will be correct for 
most macroscopic experiments. An important exception might be given by an experiment of 
nanoscale particles in glasses pQ. In that case ag and bo have the same magnitude, but Db 
seems to be one order of magnitude larger than Da- This leads to cq S> bo in contradiction 
with Eq. (|15p . However, since only one experiment was carried out yet, we cannot see whether 
Eq. |H]) holds. 



3.3.4 The width law 

The width law has been neglected in the discussion of Liesegang pattern formation for a long 
time because of difficulties in defining and measuring the width w n of the bands consistently. 
To obtain a finite width in microscopical models one needs to introduce a non-local growth 
term and a second threshold as N(c, d) from Eq. (|T3]) . On the other hand it is also possible to 
generate a width using Eq. IT2|) for macroscopic models, e. g. [2T] , 

The first theoretical step in the direction of clarifying the law was done by Droz et al. [21] 
using a model based on Eq. <|T2|) . They showed that the total mass m n of D material in the nth 
band is proportional to x n . The concentration profile d(x,t) for t ^> t n+ i will be stationary and 
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has a scaling function of the form, d n (x,t) — A n d[(x — x n )/w n ] for x n < x < x n + w n , where 
the amplitude A n could depend on n. It is thus possible to relate m n with d n by integration, 
m n = A n J* n+Wn d\^x — x n )/w n ] dx — jA n w n oc x n , where 7 is the substituted integral over 
d n . Using the first version of Eq. Q, w n ~ x", this yields A n ~ xf^ with a + f3 = 1 [24]. The 
authors also discussed experiments, fitting the results with 0.9 < a < 1. 

A more direct treatment was proposed by Racz 23 . Because the A + B — ► C process is 
independent of the C — > D process the number of C particles produced in the first process is 
the same as the number of D particles ending up in the bands. The number of C particles can 
than be calculated using the assumption that the constant concentration Cq behind the front 
is really reached. Later, C precipitates into bands of D with a high concentration dh while 
the low concentration cj of C remains between the bands. The particle conservation law yields 
)co = (x )q + w n d h leading to the width law in the form 

co-Q c , 

W n = V~, X n S3 p — X n , (17) 

dh - ci dh 

since c; is usually very small. If dh varies with the band position as dh ~ x^ the width law 
takes again the form w n ~ x" with a + = 1. 



3.4 Further reaction-diffusion models and comparison 

For charged C particles it is possible to extend the nucleation and growth model such that 
nucleation only occurs if the concentration of the outer electrolyte A exceeds a third thresh- 
old K a . Such models are known as induced sol coagulation models [11 55|56| . The additional 
threshold is motivated by the fact that the repulsive electrostatic interaction between the C 
ions can be screened by A particles. The model yields band formation significantly behind the 
reaction front as observed in some experiments [12] . The functional form of p deviates from 
Eq. O [32], 

-^(l-K a /a )] - 




(18) 



A fourth category of models is called competitive growth models [1615715 8 59 60 61162] . In 
these models, the precipitates D can dissolve with a probability decreasing with increasing 
size of the precipitates. In special cases it is possible that the bands move [53]; dissolve and 
reprecipitate [TS], such that the total number of bands leads to a chaotic time series. 

It is possible to compare the F and G functions of the Matalon-Packter law for the 
supersaturation models with thresholds discussed in Sects. I3~2l to 13.31 and Eq. (fT5)) . As already 
mentioned, the law is valid only for bo/a <C 1 and Db/Dj < 1. The ion product supersaturation 
model (see Sect. 13.21) predicts F(b ) ~ const and G(bo) ~ K sp /b^ : which is in contradiction 
to experimental results yielding monotonously decreasing functions of bo for F and G [8] . The 
nucleation and growth model predicts F(b ) ~ G(b ) ~ K sp / (ab — K sp ) w K sp /b if K sp <C bo. 
The induced sol coagulation model predicts F(bo) ~ Kl p jb\ and G(bo) ~ ( a /bo + P/bg). In 
addition, Antal et al. }32j remark that a refined version of the sol coagulation model will 
converge to the nucleation and growth model in the limit K a /ao — > 0. With such a model 
it should be possible to vary the functional dependence of F such that F(b ) ~ 6q 7 with 
1 < 7 < 2, corresponding to 0.2 < 7 < 2.7 observed in experiments, see [32] and references 
therein. Such refinement is only possible if the C particle are charged. For uncharged C particles 
the nucleation and growth model still yields the best predictions. 

Comparing all four supersaturation models one can conclude that the nucleation and growth 
model serves best as a reference model. It is the easiest model yielding the Matalon-Packter law. 
Although it needs an intermediate state C, the model is actually most suitable for analytical 
and numerical studies. The sol coagulation model and the competitive growth models, on the 
other hand, can explain details observed in specific experimental setups, but they are not needed 
to explain the basic universal laws discussed in Sect. [2] The variety of models yielding these 
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Fig. 4. (colour online) Illustration of the spinodal decomposition model. The homogeneous part of the 
free energy F is shown as a function of the 'magnetization' ra = c — (c; + dh)/2. See text for detailed 
explanation. 

laws is remarkable. The same is true for the variety of physical, chemical and even geological 
systems showing the phenomenon. This fact suggests that a very general mechanism governs 
the pattern formation. Such a mechanism will be described in the next section. 



3.5 The spinodal decomposition model with Cahn-Hilliard dynamics 

Although the supersaturation models described in the previous sections account for most of the 
Liesegang phenomena, there are drawbacks. Firstly, the threshold parameters controlling the 
growth of the bands are difficult to grasp theoretically and not easy to control experimentally. 
Secondly, it is difficult to derive how the band formation could be manipulated in a desired way, 
since the structure of the models is too complicated. Furthermore, the specific models might 
seem insufficiently universal. Hence, a new approach free of thresholds and reaction-diffusion 
equations was recently suggested by Antal et al. [27] . They studied only the second process 
C — > D of Eq. (|5|) applying the spinodal decomposition theory for phase separating processes 
64 65J to describe the phase separation into bands. We note that the distinction between C and 
D is actually not necessary here, since the bands correspond to areas with high concentration 
c = dh, while there is little C between the bands, c — ci (see also Sect. 13.3. 4[) . Although the 
model was introduced to describe C — > D processes it is equally valid for processes where the 
C particles arrange in bands of high and low concentrations. Then the phase separation is 
assumed to take place at a very low effective temperature leading to stable bands after long 
times. 

As we have seen in the Sect. !3.3Tl the production of C particles in the A + B — > C reaction 
is well understood. The reaction front - an approximately Gaussian shaped region where C is 
produced at the rate R(x,t) given by Eq. - moves diffusively with its centre at Xf(t) = 
^/2Dft and its width increasing as wj(t) ~ t 1 / 6 . Hence, the spinodal decomposition model 
can start with the C particles. Their dynamics are described by a simple phase separating 
equation taking particle conservation into account. The specific dynamics were introduced by 
Cahn and Hilliard [66167] . We note that an equivalent approach is the so-called model B in 
critical dynamics [68] . 

In the Cahn-Hilliard equation the concentration of C, c, is represented by the so-called 
'magnetization' 

m = c - (ci + d h )/2. (19) 

Then a Ginzburg-Landau-type free energy is defined in the simplest possible way required 
for obtaining two minima. This free energy, illustrated in Fig. [U is finally inserted into the 
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Cahn-Hilliard equation [23 , 

d t m = -\d 2 x [em - 7m 3 + ad%m] + R(x, t), (20) 

Here, R(x, t) is the source term introducing new C particles into the system via the A + B — > C 
reaction. The parameters e and 7 have to satisfy \J e/7 = [du — c;)/2, while a must be positive 
to eliminate short- wavelength instabilities. A and a can be used to set the time-scale and length- 
scale of the system, which leaves e as the only free parameter. Since e measures the negative 
deviation of the temperature T from the critical temperature T c , e > is needed for T < T c . 
No phase separation occurs for T > T c . 

The model can be understood most easily by looking at the free energy sketched in Fig. 2] 
The source term R(x, i) moves the system form m = —m e (for c = q w 0) over the spinodal 
point — m s to m = mo which corresponds to c = cq. Choosing mo such that it is located in the 
unstable regime, the system will move to the second stable regime at m = m e corresponding 
to the high concentration c = dh, i- e., band formation. The band becomes a sink for C, and 
its width grows until the front moves away so that m can decay below m s , and the unstable 
regime is reached again, moving the system back to m = —m e . Using this model it is possible 
to produce Liesegang patterns satisfying the spacing law ([2]) in agreement with the Matalon- 
Packter law [23127169] , Due to the conservation of C particles the arguments regarding the 
width law presented in Sect. 13.3.41 also apply here. Actually, the derivation of the width law in 
Eq. (fl~7|) was introduced in the context of the spinodal decomposition model. 

A special feature not observed in the threshold models is the low density phase with non-zero 
c, which has been reported in many experiments. The spinodal decomposition model can also 
be implemented for sophisticated conditions. For example, the effect of an additional electric 
held was discussed recently this way [42J . However, this can be done similarly using a threshold 
model (70). 



3.6 The kinetic Ising model with Glauber and Kawasaki dynamics 

An alternative way to model the phase separating dynamics was recently proposed by Magnin 
et al. [55] along the lines of the kinetic Ising model for ferromagnets. The main advantage is 
that this model fully describes the fluctuations, going beyond the mean-held approximations. 
Like in the spinodal decomposition model (see previous section), the concentration of C (and 
D) particles is represented as a magnetization. The particles are identified as spin-up sites in 
a cubic lattice, while spin-down sites represent vacancies. Then the formation of precipitates is 
modelled by a combination of spin-flip and spin-exchange dynamics |71j . The Hamiltonian is 
the usual nearest neighbour Ising Hamiltonian with ferromagnetic coupling J > between the 
spins <7 r at site r, modelling the attraction of the C particles, 

H = -J a ^'- ( 21 ) 

neighbours r,r' 

Specifically, Glauber dynamics [72| are used to add the C particles in an initial state, in 
which all spins are down. The spin-flip rate w r at site r is given by w r — R(v, i)(l — c r )/2, with 
R the source term discussed in Sect. 13.3.11 (see Eq. ©). Since the width Wf(t) ~ i 1 / 6 of the 
reaction front is not changing much, a constant width leaving behind a constant concentration 
Co of C particles was chosen |28) . The diffusion and interaction of the spin-up sites, i. e., of 
the particles, is modelled by a spin-exchange process with Kawasaki dynamics [73,, w r ^ r i = 
[1 + cxp (SE I '(ksT))]" /r e . This exchange rate satisfies a detailed balance at temperature T. 
The flip frequency r e sets the time scale, T is the temperature, ks the Boltzmann constant, 
and SE the energy change. 

2d and 3d simulations have been reported for this model 28J. In 2d no pattern formation 
could be found for a wide range of parameters, in contradiction with previous work using lattice- 
gas simulations [22l25j . In 3d patterns emerge in a restricted parameter range. Due to limited 
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computational power no quantitative results for Liesegang patterns have been published yet 
and fluctuations have not been analysed, indicating that the model needs further investigation. 

In the next section, we will thoroughly discuss lattice-gas simulations for Liesegang pat- 
tern formation, since these studies have already yielded quantitative results and fully include 
fluctuations. 



4 Lattice-gas simulations 

In the previous section, we have seen that several mean-field models can reproduce the basic 
laws of Liesegang pattern formation. However, there are problems not solved by mean-field 
models. 

1. Is the spacing law valid for all distances or only asymptotically? 

While simple mean-field theories yield x n cx (1 + p) n or equivalently x n +i = (1 +p)x n , see 
Sect. 13.21 experimental works usually report an asymptotic behaviour x n+ \ — > (1 +p)x n for 
large n (Eq. ©). 

2. Which version of the width law is the true one? 

While the mean-field theories generally yield w n oc x" (see Eq. (|17|l ). experiments have been 
fitted successfully by both version of Eq. ((4|. 

3. Which deviations from the Matalon-Packter law {3J) or its more general mean-field form 
{1J$ are relevant? 

Although most experimental findings are consistent with the Matalon-Packter law ([3]), one 
can expect deviations if Db is large or if fluctuations become important. 

4. How stable are Liesegang patterns forming under different conditions and in nanoscale sys- 
tems? 

Simple mean-field models assume the existance of quasi-periodic Liesegang patterns and thus 
cannot be used to study the stability of the pattern formation process. In general, mean-field 
quantities like concentrations are not well suited for studying nano-sized systems, since the 
number of atoms of one agent in a given small volume may fluctuate significantly. These 
fluctuations can increase or decrease the stability of the Liesegang patterns. 

5. What Liesegang patterns can be expected under special, restricted (e. g., dimensionally re- 
duced), or designed geometries? 

This question is particularly important if self-organization of Liesegang patterns shall be 
applied to design nanoscale devices. However, mean-field models usually assume a quasi- Id 
geometry. 

In this section we will address problems 1 to 3 by means of lattice-gas simulations. Problems 4 
and 5 will be studied in later publications; we just briefly comment on them here. 

Problem 4 can be addressed only by models that include the fluctuations of the particle 
numbers and thus go beyond the mean-field limit. Two major approaches in this direction have 
been published so far: the kinetic Ising spin simulation reviewed in Sect. 13.61 [28 . and lattice-gas 
simulations on the basis of the nucleation and growth model by Chopard et al. [22|25j . Since 
the Ising model approach is still not so far advanced, we chose the second approach here. An 
additional advantage is the existence of a corresponding mean-field model (see Sect. 13. 3| ). that 
the results can be compared with to find out deviations in the universal laws and effects of 
fluctuations. Equations (fT4"]) and (jTTJ) represent the mean-field solutions for the spacing law and 
width law, respectively. 

Problem 5 has already been addressed by Liesegang himself since the first experiments 
were in polar geometry. However, most theoretical work was done in quasi-ld geometry, and 
nobody investigated how geometry affects the pattern formation and the empirical laws. A first 
step towards new geometries was done in [2] using mean-field models plus a stochastic term. 
However, this ansatz seems to be hard to control for some special geometries. Therefore we 
think that the lattice-gas model can serve as a good candidate to test also new geometries and 
their influence on the pattern formation. 

Lattice-gas simulations can serve as a computational experiment to test how the mean-field 
solutions can be applied to 'experimental' data and furthermore how fluctuations might cause 



14 



deviations from these solutions. As suggested by the separation of the two processes in the 
nucleation and growth model, the lattice-gas simulation consists of two stages. In a first stage 
the C particles are generated (cf. Sect. 13.3. ip . The second stage simulates the precipitation of 
the C particles by rules comparable with cellular automata (cf. Sect. 13.3. 2j) . 



4.1 First process: reaction-diffusion A + B — > C 

In the simulation, we consider a simple cubic lattice of size M X M X L with L ^> M, see 
Fig. 2(a). Each lattice site represents one cubicle (cell) of the system. Initially, the lattice is 
homogeneously filled with B particles, i. e., there are on average &o independent B particles on 
each site. The A particles are placed on the left plane of the lattice with ao particles on each of 
the M x M sites. The parameters ao and bo can thus be interpreted as concentrations per lattice 
site. Typically, bo is in between 10 and 200, while ao is about twice as large. The mean-field 
limit can be reached if either the considered cells are enlarged or the number of particles per 
cell is increased. Thus, increasing both ao and bo but keeping their ratio constant, drives the 
simulation to the mean-field limit. 

The dynamics is modelled as follows. To simulate diffusion, both, A and B particles per- 
form independent random walks on the lattice. The diffusities Da and Db are defined as the 
probabilities that a motion in either of the six possible directions takes place in a given time 
step |74|75j . i. e., Da and Db are proportional to the physical diffusivities. After each time 
step, the left plane of the lattice is re-filled with A particles with the initial condition. A re- 
action A + B — > C takes place with probability k (k = 1 in our simulations), if at least one 
A and one B particle is found on the same site |51|76j . Afterwards, all C particles also diffuse 
independently on the lattice with diffusity Dq- It was previously shown that the mean-field 
predictions regarding the reaction rate R(x,t) (reviewed in Sect. [3~3~Tj) are in agreement with 
the simulations in 3d [51] . while logarithmic corrections occur in 2d [76] . 

Here, we use for the first time 3d lattice-gas simulations for Liesegang patterns in contrast 
to the 2d setup employed by Chopard et al. [22]. For an alternative 3d simulation see Sect. 13.61 
|28j . To confirm that the mean-field limit is reached for large concentrations, we ran simulations 
with different bo (keeping ao/&o constant) and found that co/bo becomes constant in agreement 
with Eqs. (fTU]) and (fTT]) . Figure [5j^a) shows the simulated deviations from the mean-field limit, 
which decay below one percent for bo 3> 2. We focus on bo > 10, where the deviations are below 
10~ 4 . In a 2d (or Id) simulation, fluctuations would alter cq such that C0/&0 does not converge 
to the mean- field limit even for large bo and ao [76] . 

Since the mean-field solutions are valid for a 3d setup, it is possible to insert the C particles 
directly. This ansatz was first used in a kinetic Ising scenario [28] (see Sect. I3.6[ ). We employ this 
idea to speed up our numerical calculations, since fully simulating the reaction-diffusion process 
A + B — > C reduces the computational speed by at least one order of magnitude. Otherwise 
we could not work with concentrations bo > 10 in sensible time. However, this approximation 
reduces the fluctuations of C particle production, in particular for small concentrations. 

The probability to insert C particles into the lattice is given by the reaction rate R(x, t) 
(see Eq. (|9])), with a centre position Xf(t) moving as x/(t) = y / 2Dft to the right. Usually, 
Df is calculated via Eq. ([It?]) , assuming that the A particles are distributed homogeneously for 
x — > — oo, i. e., a(—oo,t) = a_oo. In this case the A particle concentration outside the quasi- 
stationary region (see Sect. l3.3TTj) can be approximated by a(x, t) — a_oo— a/ 

with a constant a/. However, in our configuration, Eq. ([T7j|) needs to be modified, since the 
concentration of A particles is held constant for all times on the left plane, i. e., for x = 0: 
a(0,i) = a-co — a/ = ao- A straightforward derivation similar to the derivation of Eq. (|10|) in 
[52] leads to the solution 

erf exp = 

\V2D A J P \2D A ) boVD^ 




2D B 



(22) 
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Fig. 5. (colour online) (a) Deviations of the 3d lattice-gas simulations from the mean-field limit. The 
values of simulated ratios Co°/&o have been calculated from the plateau of c(x, t) observed behind the 
diffusion front (see Fig.[3jb)). The corresponding mean-field value Cq F /6o was calculated from Eq. {TT}. 
The red curve marks the (scaled) standard deviations of the particle numbers per lattice cell from their 
means Cq G indicating the fluctuations in the production of the C particles. The simulation parameters 
are ao/&o = 2, Da = 1, Db = 0.1. (b) Results of our lattice-gas simulations with C particles inserted in 
an approximated reaction zone according to Eqs. (|23p and (|24p . The values of simulated ratios Cq G /co 
have been calculated as in (a). They are in good agreement with the inserted Co for all concentrations. 
Red curve same as for (a). The parameters are Df = 1.22 (calculated from Eq. ([22}) and A — 3. 



Following [28] we approximate the nearly Gaussian-shaped reaction rate term R{x 1 1) by 

R(x,t) = -j=e(x-Xf + A)0(xf +A — x), (23) 
yt 

with constant width 2A of the front. These simplifications are justified since the band forming 
process does not depend on the exact form of the reaction zone 27j and the width of the 
front increases very slowly in time (as i 1 / 6 , see Eq. ©). The most important feature of the 
reaction front is that it leaves behind a constant concentration Co of C particles. Therefore, Af 
in Eq. (f2"3"| is chosen to be 

Af = \^c . (24) 

Figure [5](b) compares the simulated Cq with the Co value inserted into the simulation via 
Eq. (f24l) . In contrast to a full simulation of the A + B — > C reaction (see Fig. E{a)), Cq G 
reaches the mean-field value already for small concentrations since the fluctuations induced by 
the motion of the A and B particles are eliminated. Only the standard deviation (red curve) of 
the number of particles per cell is similar as in Fig. [5Ja), since the reaction front in (b) is an 
approximation of the reaction front in (a) and both cause the same fluctuations in C particle 
production. 



4.2 Second process: precipitation C — > D 

The second part of our lattice-gas simulation is the precipitation of C particles, generating the 
immobile precipitate D. The nucleation of D depends only on c(x,t), while the growth of D 
depends on both, c(x,t) and d(x,t), see also Sect. 13.3721 In the particle picture of our lattice-gas 
simulation the density will be the number of particles divided by the considered cell volume. 
Here, this considered volume around a given lattice cell includes all neighbour cells, i. e., the 
27 cells in a cube of 3 x 3 x 3 cells. This definition sets the mean reaction distance. 

Two thresholds are introduced. If the mean local concentration of C particles (in the 27 
cells) exceeds a threshold K sp , nucleation occurs. C particles in the vicinity of D particles 
precipitate already if their mean local concentration exceeds a threshold K p < K sp . C particles 
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Fig. 6. Spacing law for simulated data depicted in a way different from Fig. [2jb), where x' n+1 /x' n was 
plotted versus n to observe the asymptotical behaviour (cf. Eq. (|26p ). Here, x' n+1 is plotted versus x' n 
to observe the linear behaviour with slope 1 + p and offset p£ (cf. Eq. (|27| l1. yielding p — 0.066 and 
£ = 62. Parameter set for the 3d lattice-gas simulation: ao = 130, b = 65, Da = 1, Db = 0.1 (leading 
to Df = 1.22 and c /6 = 1.072); D c = 0.1, K sp /b = 0.93, and K p /b = 0.52. 

on top of D particles always precipitate, which makes the growth process fast enough to deplete 
a region of C particles. Without this option the growth process would not terminate, and no 
additional Liesegang bands could be formed. The resulting bands have distinct positions, which 
we will denote by x' n in the following; their width is denoted by w n . We calculated x' n as the 
centre of the bands, i. e., the mean of the positions of the first and last D particles within each 
band. The corresponding widths w n are defined as the differences between these two positions. 

4.3 Reconsideration of the spacing law 

Equation @, i. e., the asymptotical convergence of the ratio of the positions of neighbouring 
bands x n+ i/x n — ► 1 + p for large n, was estimated from empirical findings; see Sect. [2j On 
the other hand, theoretical analyses suggest that the ratio x n+ i/x n is identical with 1+p 
[29I30I31I46) . as was discussed in Sects. 13.21 and 13.31 Although these two forms of the spacing 
law contradict each other, both are still used in parallel without much considerations. We 
propose a way to reconcile the experimental findings described by the empirical form of the 
spacing law with the theoretical mean-field prediction. 

In experimental and numerical results the position of the first Liesegang band is usually 
blurred and thus not well defined, contrary to theoretical analyses. Thus, there is an arbitrari- 
ness in choosing the position of the first band. Therefore, it is better to use different variables 
x' n for the experimentally or numerically measured band positions and x n for the theoretical 
positions, since there may be an offset £ between them, 

x n = x'„ + £, (25) 

due to, e. g., the blurred first band. Another possibility is that the mass of D in the first 
Liesegang band as well as the band's width cannot be close to zero as would be required if the 
linear increase observed in Fig. [^c) was beginning at x' = 0. We will see later that usually 
£ > 0, i. e., x n > x' n , indicating that the ideal laws are based on an 'imaginary' starting point 
outside the real sample. In literature, however, significant confusion is caused by the fact that 
both variables, x n and x' n , are not distinguished. 

Studying measured values x' nl and assuming both, x' n — x n — £ and the ideal mean-field 
result, x n+ \ = (1 +p)x n , one finds the asymptotic (empirical) form of the spacing law, 
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Fig. 7. (colour online) Reconsideration of the Matalon-Packter law based on 3d lattice-gas simulations, 
(a) The yellow squares indicate simulation results for the spacing law coefficient p obtained keeping bo 
constant and varying ao- The straight line is a linear fit to the first four data points using Eq. (|3]). (b) 
Simulation results similar to those shown in (a). Df and Co are calculated using Eqs. (|22[) and (|ll|l . 
respectively. The blue curve is a fit of Eq. (|14|. while the red curve takes fluctuations into account. 
Parameters: b = 55, D A = 1, D B = 0.1, D c = 0.15, K sp /b = 0.96, and K p /b = 0.52. 



for large n as in Eq. ((2]). Figure Efb) shows this asymptotic behaviour as observed in simulated 
data; note that actually the values of x' n rather than x n are plotted. The coefficient p is hard 
to determine in such plots, especially if only few bands are present. A more convenient way to 
extract p from measured data is to fit the equation 

x' n+1 = (l+p)x' n +pt, (27) 

where p£ is the intersection with the ordinate, as shown in Fig. [5] for the same data. Even 
for a small number of bands, p and £ can be extracted conveniently. In addition, we believe 
that the offset £ is an important ingredient that should not be neglected in interpreting any 
experimental or simulated data. However, it remains unclear if or how the value of £ could be 
predicted. 



4.4 Reconsideration of the Matalon-Packter law 

The original Matalon-Packter law ([3]) predicts a linear dependence between bg /ao and P- It was 
shown in Sect. [373731 that this empirical law holds for &o <C ao only. In addition, a more general 
law was presented in Eq. (fT4"]) [32] . The deviations from Eq. ([3J are even stronger in the results 
of lattice-gas simulations. An example is depicted in Fig. [7][a). A linear dependence is observed 
only for ao/fro > 250, a criterion shown to be equivalent to Db/Dj « 1 [32]. If ao is reduced for 
constant bo, Df increases, driving the system out of the regime in which Eq. ([3]) is valid 32J. 
Our lattice-gas simulations confirms these results, see Fig.[7][a). 

Figure [7][b) compares a fit of the generalized Matalon-Packter law Eq. (TH)) (blue curve) 
with the results of our lattice-gas simulations. The deviations are still quite large. Only if 
fluctuations are taken into account by further modifying the generalized Matalon-Packter law, 
a nice agreement can be reached. The formula used for the red fit in Fig. [T^b) will be discussed 
and motivated in detail in a later publication. 



4.5 Reconsideration of the width law 



As discussed in Sect. [2] two competing empirical forms of the width law are used for fitting 
experimental data, 

w n = ii\x n + /i 2 and m„ a (jl]) 
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Fig. 8. Reconsideration of the width law based on 3d lattice-gas simulations, (a) Linear presentation 
of w n versus x' n (points) together with a fit of the first version of Eq. ((4]) with /ii = 0.035 and fi2 = 3.4. 
(b) Double logarithmic presentation of the same data together with a fit of the second version of Eq. Q 
with slope q = 0.82. The same simulation as for Figs. [2] and [6] was used. 



Figure [5] shows the results of our lattice-gas simulations in both representations, because w n is 
plotted versus x' n both linearly and double logarithmically. Clearly, it is not possible to favour 
either of the two forms of the width law, since all widths are small compared with the size of 
the sample and both fits have an equivalent quality. 

On the other hand, in Sect. I3.3T41 a simple and general width law was deduced theoretically 
(see Eq. ([T7|) ) using particle conservation. Assuming scaling behaviour of the band densities d n , 
d n oc Xfl, and setting q = (since no precipitate occurs between the bands in our simulation), 
one obtains 

w n = p-j-x n OC x% (28) 

with a = I — p. This law was derived using the theoretical x n where the ratio for two consecutive 
bands is constant, i. e., the exact form of the spacing law. A width law valid for experimental 
or numerical data of the positions x' n of the bands can be derived introducing the offset £, see 
Eq. (|23|) . The analytical form of the width law thus becomes 

Wn = V-T- {x' n +0= P-T- X n + (29) 

which is in between the two competing empirical forms Q, since d n oc {x' n + O - Note that 
p£ is identical with the offset of the fit for the spacing law (|27[) . Hence, the parameter £ is also 
important for understanding the width law. 

To test these predictions we have run simulations and fitted Eq. ()28|) to extract a from the 
data disregarding £, i. e., inserting x' n directly for x n - The results are shown in Fig. [9]^ a), yellow 
points. In a second approach we used the same data and inserted x n — x' n +£ into Eq. (|28|) . also 
obtaining a as a fit parameter. The results are shown in Fig.[9ja), red points. Both approaches 
lead to different values of a indicating that the offset £ plays a significant role. 

Since the width law depends on the density of the bands d n , it is not possible to distinguish 
the correct form by just looking at the width exponent a. Therefore, we have also calculated 
the densities d n in the numerical simulations. Figure [UJb) shows the exponent f3 obtained by 
fitting d n oc x' n (yellow points) and d n oc (x' n +£)' 3 (red points). In addition, Fig.[9lc) depicts 
the sum of both exponents from parts (a) and (b), a + f3. The product w n d n can be interpreted 
as the mass m n of band n; it scales as m n = pcox n — pcox' n +pco£ according to Eq. (|29|) . Thus, 
according to theory, the mass exponent a + (3 must be one. The results shown in Fig. [9^c) 
indicate that this holds only in the case where the offset £ is taken into account correctly (red 
points), confirming that £ must not be disregarded. 



Comparison of models and lattice-gas simulations for Liesegang patterns 
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Fig. 9. (colour online) Fits of (a) width, (b) density, and (c) mass scaling behaviour based on our 
3d lattice-gas simulations for different values of Co. The exponents have been obtained by fitting a 
power-law to plots of (a) w n , (b) d n , and (c) m n = d n w n versus x' n (yellow dots) and versus x' n + f 
(red dots). The parameters of the simulations are identical with those used for Figs. [2] and [6] except 
for varied co. 



4.6 Application to previous experimental data 

We believe that much confusion in literature could be avoided if the difference between x n and 
x' n was taken into account. In particular, the long-standing discussion about the width law 
could probably be solved if the analysis of experimental data included plots versus x n = x' n + £- 
We suggest that former experiments should be re-analysed using £ and x' n to confirm our 
conclusions. 

An experiment which confirms the result of Eq. (129)) was published recently [19] . The authors 
investigated pattern formation in a K-Carrageenan gel. They measured x' n and w n for different 
concentration of the outer solution. First they plotted x' n+1 versus x' n in a fashion similar to 
Fig. [6] and extracted p. Looking at the plot, we can clearly see offsets corresponding to p£. 
Secondly they plotted w n versus x' n in a linear fashion similar to Fig. [5Ja), extracting the 
parameter /ii = pco/d n in Eq. (|2Uj). assuming constant densities d n . Looking at the plot, we 
also clearly see offsets corresponding to fi 2 = p£,c /d n (according to Eq. (I29)) ). In addition they 
determined co/d n = 0.53. This means that the offsets of the width law must be approximately 
half of the offsets of the spacing law. This conclusion is in quantitative agreement with the 
offsets we read from the plots in |19) . The experimental data thus confirm our conclusion that 
the offset £ should be taken into account. 

The differences of Eqs. (fl~7|) and (f29|) explain the different approaches used in experiment 
and theory. In summary, we have solved problems 1 and 2 raised in the beginning of Sect. |4j 
problem 3 was discussed in Sect. 14.41 



5 Summary and outlook 



In this paper we reviewed the four empirical laws believed to govern Liesegang pattern formation 
and several important mean-field models reproducing these laws as well as a few Monte-Carlo- 
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type simulations. Based on our detailed 3d lattice-gas simulations of the nucleation and growth 
model, wc detected three major problems in reconciling the experimental reports with mean- 
field results. However, our simulations also helped us to find a straightforward solution. 

We have shown that all basic empirical laws describing Liesegang pattern formation, i. e., the 
time law, the spacing law, the Matalon-Packter law, and the width law can be understood on the 
basis of the nucleation and growth model. A full agreement between experimental observations, 
simulation results, and mean-field models can be obtained only if a constant offset £ between 
measured and theoretically assumed band positions is taken into account, as suggested in this 
work. A re-analysis of a previous experiment concerning the spacing law and the width law 
confirmed our suggestion. We propose that further experiments should be re-analysed to test 
the refined predictions. 

In addition, we suggest further experiments with systematically varied concentrations of 
both, inner and outer electrolyte to confirm or extend the generalized Matalon-Packter law (fl"4|) . 
Our current work in progress regarding the effects of fluctuations on Liesegang pattern formation 
indicates that additional modifications of the Matalon-Packter law are necessary in small-scale 
systems. Reliable experimental results for the nanoscale regime, where fluctuations are definitely 
important, are expected to become available soon, since experiments with nanoscale particles 
get feasible. A big open question will be how to manipulate the pattern formation such that 
interesting nanoscale devices could be designed. A deeper understanding of the dynamics is 
important for reaching this goal. For a first promising work in this direction, see [43j . 

We thank I. Bena, M. Droz, M. Dubiel, I. L'Heureux, Y. Kaganovskii, I. Lagzi, K. Martens, and M. 
Rosenbluh for discussions. We would like to acknowledge support of this work from the Deutsche 
Forschungsgemeinschaft (DFG, project B16 in SFB 418). 
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